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Abstract 

Belle data for the tttt mass spectrum in X(3872) — > tt + tt^J/^ 
are well fitted by p + J/^ with J PC = 1 ++ , but are poorly fitted 
F ^ , by the tttt S-wave. This rules out quantum numbers 1 H , 2 and 

pL^ ! 3 . Formulae for partial wave amplitudes are given for all likely 

• J PC and decay modes of the 3941 MeV peak which Belle observe in 

CN ■ to + J/^>. Angular correlations involving up to five angles may allow 

\ a spin-parity determination, even with quite low statistics. Special 

1 attention is given to the case where X(3872) is a cusp or quasi-bound 

state with the same quantum numbers as the 3941 MeV peak. 



PACS Categories: 11.80.Et, 13.20.Fe, 13.20.He, 14.40.Lb 



1 Introduction 

Oh 



The Belle collaboration observes a narrow peak at 3872 MeV in tt + tt~J/^; 
it has a width < 2.3 MeV with 95% confidence [1]. This peak is confirmed 
^Cj . by CDF II [2] and by the DO [3] and Babar [4] collaborations. There has been 

extensive discussion of possible quantum numbers and expected charmonium 
levels [5-14]. The fact that X(3872) lies very close to the D°D*° threshold 
and is narrow has prompted suggestions that it may be a quasi-bound state of 
DD* + D*D [15]. Tornqvist ascribes this possibility to long range attraction 
due to ii exchange [16]. A possibility explored here is that it may be a cusp 
due to de-excitation of DD* to other open channels [17]. 

A new result is that Belle also observes an enhancement near threshold in 
luJ/^ at 3941 ±11 MeV with a width r = 92 ±24 MeV [18]. As a shorthand, 
this peak will be called K(3941). 

This paper discusses three somewhat disparate topics related to quantum 
numbers J PC of X(3872) and F(3941). Section 2 examines fits to the tttt mass 
spectrum in X(3872) — > 7t + tt~ J for all likely values of the orbital angular 
momentum L between J/\I/ and tttt. For L = 0, the p(770) gives an excellent 
fit, as Belle remarked in their first publication [1]; this supports J PC = 1 ++ . 
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Fits with the tttt S-wave require discussion of how to parametrise it; even the 
most conservative parametrisation gives distinctly poorer fits because S-wave 
structure is broader than p(770). Discrimination against quantum numbers 
1+- O-Pi), 2— ( 3 D 2 ) and 3~ ( 3 D 3 ) is quite strong. 

Section 3 reviews remaining possibilities. The V(3941) lies close to pre- 
dicted radial excitations of XcOi Xci an d Xc2- However, the observed decay 
width of 92 ± 24 MeV is a factor ~ 10 larger than expected for uoJ alone. 
The remaining width is likely to come from decays to DD* or DD. Predic- 
tions of Eichten et al. [14] then favour J PC = 1 ++ . One is then faced with 
the possibility that X(3872) and F(3941) may both have J PC = 1++. This 
can arise if X(3872) is a cusp at the DD* threshold or a 'molecular' state. 
In the latter case, mixing between 1 ++ S and D-waves needs experimental 
study. 

Progress depends on angular momentum analyses of both X(3872) and 
F(3941). Section 4 is intended to assist such analyses by providing tensor 
formulae for partial wave amplitudes. Up to now, discussions have depended 
on decays of X or Y. However, there is additional information available 
from the production process B — > K + X (or 7). The full formulae involve 
correlations between five angles describing (i) this production process, (ii) 
decay angles of X (or Y) — > C + D and (iii) decay angles of C and D. 
Tensor expressions have the virtue of being written compactly in Cartesian 
coordinates and can be programmed in a few lines of code. These formulae 
extend naturally to DD* decays and indeed most decays of the B meson 
and charmonium states; only minor substitutions of variables are required. 
Hopefully these formulae may find application to a large range of processes. 
An Appendix discusses the three-dimensional geometry of the J/\I/ decay and 
the effect of the Lorentz transformations to the rest frames of X or Y. 

Earlier, Pakvasa and Suzuki have given angular distributions for some 
decays of X(3872) [19]. Rosner [20] gives decay angular distributions for 
important cases: J p = + , 1 + and 2 + S-wave decays to (p° or a;) + J 
and for J p = 0~; he does not discuss correlations between decays and the 
production process. Full formulae for both production and decay are given 
here for all J p up to 3~. 

Section 5 summarises conclusions and makes some suggestions about the 
way angular momentum analyses might proceed. 
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2 The 7r + 7r mass spectrum of X(3872) decays 

Belle point out that this mass spectrum peaks at the highest available mass, 
close to p(770). A second possibility which needs quantitative discussion is 
interpretation in terms of the mr S-wave. This is a tricky point. 

Fig. 1(a) shows the Breit-Wigner line-shape for p(770) and the intensity 
for 7T7T S-wave elastic scattering, denoted by S. Both intensities have peaks 
close to 770 MeV, but the S-wave intensity is much broader. This figure does 
not include the phase space available for X — > nirJ/^. 
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Figure 1: (a) Line-shapes of p(770), nn elastic scattering S, and the a pole 
a. (b) The nil mass spectrum from Belle data on X(3872) — > 7r + 7r~J/\I/, 
compared with the expected p mass distribution (full curve), and decays to 
S with L = 0, 1 and 2. 

In BES data for J/^> — > Lun + ir~, a very different irir peak is observed at 
~ 500 MeV and is fitted with a a pole at 541 MeV [21]; the fitted intensity 
is illustrated by the chain curve labelled a in Fig. 1(a). Data of E791 on 
D + — > 7r + (7r + 7r~) are similar but give a somewhat lower mass, 478 MeV and 
a narrower peak [22]; for present purposes, the difference between these two 
is not important. What needs discussion is the difference between the a and 
S peaks. What should be expected in X(3872) -> [ttttJ/^]? 
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Actually both a and S curves may be fitted by a single formula, given 
in Ref. [23] (with a correction of a sign error in the Erratum). The reaction 
J/ty — > uiyti is a 'hard' process producing the tttt pair with a large momentum 
transfer. The Breit-Wigner denominator alone creates the peak. However, 
7T7T elastic scattering is a 'soft' process and the numerator of the amplitude 
contains an additional zero at s = 0.5m^, just below threshold. This is the 
Adler zero of Chiral Perturbation Theory. In elastic scattering, the side of 
the pole close to s = is sheared away by the Adler zero. For elastic -kit 
scattering, the intensity is then rather featureless at low mass and peaks 
only at 800 MeV. For complex s, the singularity resembles a mountain ridge 
running at 45° to the real s axis. The peak moves dramatically with Ims; 
the resonance is at 925 MeV for real s, but at 540 MeV for Ims = 0.25 GeV. 

The decay X(3872) — > SJ/^ involves a small momentum transfer to the 
7T7T system, so the nn elastic scattering amplitude S is likely to be more 
appropriate. However, it is not yet known how rapidly the Adler zero disap- 
pears with increasing momentum transfer. The approach here is to fit using 
the 7T7T elastic scattering amplitude; this is the most conservative approach, 
closest to data on X(3872). Any effect of the disappearing Adler zero will 
shift the peak downwards in mass, away from the p, making fits to the nn 
S-wave even less likely. 

Fig. 1(b) shows the 7T7T mass spectrum from Fig. 3 of Ref. [1] in 50 MeV 
bins. The full curve shows the expected signal for L = decays, folding in 
the available phase space for X(3872) — > pJ/^ . This gives an excellent fit 
with x 2 — 4.6 for 6 degrees of freedom. The mass spectrum for [pJ/^] L= i 
gives x 2 = 12.1 and cannot be excluded; however, a strong rate would be 
surprising, since the p is so close to the top of phase space. 

The dashed, dotted and chain curves show fits to Belle data using curve 
S of Fig. 1(a) and L = 0, 1 or 2. The L = fit is affected much more by 
the phase space cut-off than the p fit. It gives a poor x 2 of 19.4; J PC = l~ 
is highly unlikely if X(3872) is a charmonium state, which should appear in 
e + e~ interactions. Seth [24] points out that X(3872) could be a J PC = 1~ 
glueball; a pure glueball would not couple to e + e~. 

For L — 1 and 2, the fits are very poor with x 2 = 49 and 343 respectively. 
Blatt-Weisskopf centrifugal barriers are included with a radius of 0.8 fm. The 
fits are much worse if the a pole of Fig. 1 (a) is used with no Adler zero. For 
[SJ/^]l = i, isospin conserving decays are + ~ (exotic), l + ~ or 2 + ~ (exotic). 
For J PC = l +_ , Belle already find poor agreement with the decay angular 
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Table 1: Possible J PC for F(3941) -> wJ/>, X(3872) -> pj/^ (C = +1) or 
SJ/^ (C = —1); also possible DD* and DD decays. Those in parentheses 
are less likely. 



distribution [25]. Also Babar find no evidence for X(3872) — > r/J/^/ [26]; this 
makes 1 + - doubly unlikely. For L = 2, the J PC = 2" ( 3 L> 2 ) and 3— ( 3 L> 3 ) 
possibilities are definitely excluded. 

Babar find no evidence for charged X(3872) [27], so pJ/^f decays are 
probably isospin violating. For L = 0, all of J PC = ++ , 1 ++ and 2 ++ are 
clear candidates. The possibilities for L = 1 decays are 
(exotic), 2~ + and 3~ + (exotic); although less likely, these are not excluded 
at present. However, the rf" is expected to lie above 4000 MeV. 

3 General considerations 

Table 1 lists possible J PC for each decay; those for C — +1 are candidates for 
both X(3872) and F(3941), while those for C = -1 refer only to X(3872) -> 
SJ/ty. The notation is that s denotes combined spins of J/ty with uj (or 
p); L is the orbital angular momentum in decay channels. L = 3 decays are 
omitted as unlikely for X and Y and to CD*. 
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Eichten et al. [14] survey charmonium states likely from 3800 to 4000 
MeV and calculate decay widths. The lowest 1 D 2 , 3 D 2 and 3 D 3 are expected 
somewhat below X(3872). The radial excitations of 3 P , 3 Pi and 3 P 2 are 
predicted in the general mass range of V(3941) and are serious candidates. 

I wish to draw attention to the possibility that X(3872) and Y (3941) may 
both have J PC = 1 ++ . There are two tentative pointers in this direction for 
V(3941), though experimental proof is needed from partial wave analysis. 

3.1 The width of F(3941) 

The observed total width of V(3941) is 92 ± 24 MeV. This width cannot 
plausibly be explained by the ooJ/^ channel alone. If F(3941) is a charmo- 
nium state, the production of a decay ui requires formation of an additional 
uu + dd pair via a process involving at least two gluons. This should lead to a 
width of the same order as Xco, Xa and Xc2, i-e. 1 to 16 MeV. The remaining 
width is likely to come from DD* or DD decays. 

Eichten et al. [14] predict for the 3 Pi radial excitation a DD* width of 
150 MeV at a mass of 3968 MeV; if the width is proportional to DD* phase 
space, it scales to 127 MeV at 3941 MeV, not too far from the observed 
value. For J PC = ++ or 2 ++ , they calculate much smaller widths for DD 
decays, although there is some sensitivity in these calculations to radial wave 
functions. For J PC = ++ , they find a 40 MeV width to D S D S at 3968 MeV, 
but the threshold is at 3938 MeV, so this channel is unlikely to account for 
the observed width of V(3941). 

The agreement of the V(3941) width with the calculation of Eichten et 
al. for J PC = 1 ++ is a mild pointer. It is obviously important to search for 
decays to DD*, DD and D S D S . 

3.2 The production process 

In the Belle and Babar experiments, the V(3941) and X(3872) are produced 
via the (V — A) current. Allowed transitions are to 1 + and 1~, and to + 
via the divergence of the axial current. These are just the channels which 
have been observed in B decays. A feature of Belle data is that there is a 
strong i^x cl (3510) signal in Fig. 3 of Ref. [25]. There is also a conspicuous 
Ki(j(3770) signal [28]. Both are reached by allowed V or A transitions. The 
final state Kxco has been observed [29] and also KD S (2317) [30]. 
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However, it is difficult to estimate the strengths of forbidden transitions. 
The W boson is absorbed into the final state and the kaon is radiated with 
one unit of orbital angular momentum (in addition to orbital angular momen- 
tum involved within forbidden V or A interactions themselves). The weak 
interaction is pointlike; angular momentum transfers depend on differences 
in the radial wave functions of B and charmonium states. 

The strength of the observed X(3872) and V(3941) signals will be taken 
here as a hint that they may both have J PC = 1 ++ , though ++ is also a 
serious possibility; 2 ++ is first forbidden. 

3.3 A bound state or cusp 

If the X(3872) is a bound state of DD*, it should appear at low masses in 
that channel, in the same way as the threshold NN 3 <Si amplitude is strong 
and can be related to the deuteron pole. If the effective range is close to 
the pion Compton wave-length, the phase shift will drop through 90° at ~ 4 
MeV above the DD* threshold. One then expects a strong threshold DD* 
signal over a mass range of 5-10 MeV. 

A possibility is that the X(3872) is a threshold cusp. The way a cusp 
works is explained in Ref. [17], but will be repeated here briefly. Suppose 
the low mass DD* system de-excites to open channels; obvious possibilities 
are [St) c ]l=i, Sxco and the discovery channel rniJ/^. Of these, the first is 
likely to be the strongest, since the available momentum is well above the 
L = 1 centrifugal barrier. 

Close to threshold, the cross section for DD* de-excitation to these open 
channels follows the familiar 1/v behaviour of thermal neutron scattering. 
The imaginary part of the DD* elastic scattering amplitude fs is given by 
/c/(47r) times the total cross section. The factor k for centre of mass mo- 
mentum cancels with the 1/v factor above threshold, producing a constant 
amplitude to first approximation. However, there is a step in a at threshold 
and hence in Imfg. The real part of the amplitude is given by a dispersion 
relation: 



The step in Im fs gives rise to a peak in Re fs at threshold. Fig. 2 illustrates 
these effects from a model calculation whose parameters will be explained 
shortly. The result is a narrow peak in the intensity, Fig. 2(b). 




(1) 
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Figure 2: Schematic illustration 1 of ^possible cusp at the M ^&^ threshold: 
(a) Real and imaginary parts of the elastic D°D*° amplitude, (b) the corre- 
sponding intensity. 

This cusp appears not only in DD* elastic scattering, but in any channel 
coupled to the DD* S-wave. The scattering amplitude is an analytic function 
of s with coupling constants gi to individual channels %. A familiar example 
is a Breit-Wigner resonance for coupling of two channels 1 and 2: 



Even if there is no resonance, the s-dependence of /12 is common to all 
channels coupled to the DD* threshold and the cusp appears in all of them. 

At present, the strengths of couplings to open channls is not known. 
There is competition between decays of D* to Dir and Dj and de-excitation 
to open channels. There is also a question of whether long-range meson 
exchanges are attractive or repulsive; repulsion would shield the interaction 
close to threshold and eliminate the cusp. 

The illustrative calculation of Fig. 2 is made along the lines of Ref. [17] 
using a scattering length approximation k cot S — 1/a, where a is complex. If 
the S-wave amplitude is written as f s = {e 2l& — l)/2ik, the S-wave amplitude 
near threshold is 



Ms) 



9i92 



M 2 - s - iM[gl Pl (s) + gip 2 (s) + ...]' 



fs = 



a a + i\a\ 2 k 



(2) 



l-iak 1 + 2A;Ima + k 2 |a| 2 ' 



8 



Ref. [17] makes comparisons with pp data, where the fitted scattering 
length is Ima = 1.8 fm. For DD*, Ima will be smaller, since open channels 
are much weaker. If the F(3941) is a simple Breit-Wigner 1 ++ resonance 
with a width proportional to DD* phase space, the scattering length is 0.4 
fm and the effective range —1.5 fm. However, the imaginary part of the 
scattering length is a matter of guesswork; it need not be related to the 
F(3941) resonance. Fig. 2 illustrates the cusp in D°D*° elastic scattering 
for a scattering length taken to be 0.4 + iOA fm. The step in the real part 
represents an effective attraction, which helps bind any 'molecular' state 
such as that suggested by Braaten and Tornqvist. If the net effect of this 
attraction and meson exchanges is strong enough, a bound state appears. 

Experimentally, the X(3872) peak is at the D°D* threshold, not the 
threshold for charged particles. This charge dependence can arise from a 
combination of n, p and A\ exchanges. A large isospin violation is then 
inevitable, since the DD* phase shift changes by 90° over ~ 4 MeV; this is 
only half the mass difference between the neutral and charged thresholds. 
Isospin violation is consistent with the observed decays to pJ/^f. For this 
reason, the second cusp at the D~D* + threshold is not shown on Fig. 2, 
though it could in principle be present. The height of the cusp depends on 
the actual scattering length; the cusp gets wider as Im a gets smaller. 

A scenario worth consideration is that there is a regular charmonium 1 ++ 

state at 3941 MeV, and in addition a threshold cusp or 'molecular' state at 

the DD* threshold. If ir exchange plays a strong role, as Tornqvist suggests 

[16], it will produce a tensor interaction which mixes S and D- waves, as in 

the deuteron. In the next Section, amplitudes will be given for the D-wave 

possibility. That mixing could well be absent for a pure cusp, and might 

allow a distinction between a cusp and a bound state. An alternative is 

that the 1 ++ state is pulled down to the DD* threshold and F(3941) has 
JPC = ++ or 2 ++. 

4 Angular Momentum Algebra 

Amplitudes will be written for B — > K + V(3941), abbreviated henceforth 
as Y; those for Y -> DD* and X(3872) -> SJ/^ may be obtained by 
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simple interchanges of symbols. These expressions are easily carried over to 
a wide variety of decays of the B meson and charmonium for s and L < 
2, J up to 3. Tensor expressions will be used, following the methods of 
Zemach [31]. However, it is convenient to simplify the algebra from the 
outset by working in the rest frame of Y. In this system, orbital angular 
momentum I in the production process is described by a 3-vector for the 
kaon momentum. The orbital angular momentum between J/ty and p (or a;) 
is likewise constructed from the 3-vector J M for the J/^f momentum. This 
greatly simplifies the algebra. 

While this paper was being written, Rosner [20] produced a paper giving 
expressions for angular dependence of decays. He does not include correla- 
tions with the production process, which give additional information included 
in the tensor amplitudes written here. His expressions for decay distributions 
have been checked against formulae given here. As an aid to experimentalists, 
amplitudes will be related to the axes he chooses. 

4.1 Choice of axes 

Fig. 3 sketches vectors and angles for the case Rosner considers: X — > pJ/^> . 
A key point is that polarisation e of the J/\& is orthogonal to the axis of 
decays to e + e~ (or p + p~) in its rest frame. The leptons are shown as they 
appear in this frame. This lepton axis is adopted as my z-axis, but amplitudes 
are written in the rest frame of J/^f and p (or a;). The choice of the xy plane 
around the z direction is not important. 

Rosner chooses his X axis along the p direction in Fig. 3 and defines 
the XY plane to contain the decay pions. His Z-axis is then normal to the 
pJ/ty axis. He describes the e + e _ axis by angles 9, 0. The Appendix gives the 
components of J/\& polarisation in his axes, eqn. (26). It also gives the effect 
of the Lorentz boost from the J/^ rest frame to that of Y or X, eqn. (27). 
This result for e is needed in equations below. The Lorentz transformation 
has a small effect, so e remains approximately transverse in that system. 

The spin 1 of the p is described by the 4-vector 

M 2 — M 2 

Q = fcfa) - fcfo) - 1 2 [fc(7n) + fc(7T 2 )]. (4) 

P 

Ignoring the small mass difference between charged and neutral pions, this 
becomes the 4-vector [hi — k 2 \ fl . The time-component E\ — E 2 is small but 



10 




Figure 3: Illustration of kinematics; 9 is the polar angle of the e + and 4> its 
azimuthal angle in XYZ axes. 

easily included into the calculations. Expressions for uj decay are likewise 
described by a 4- vector W^. Amplitudes will be written here for ujJ/^f decays; 
those for pJ/^ decays are obtained by substituting Q for W. 

4.2 u decay 

The uj decay may be viewed as uj — > [/07t]l=i. For uj — > p + n~: 

W„ = e mM {K ~ h) a k p _{k + + k.+ k y = 2t mM k%k t5 _kl (5) 

the latter follows from energy-momentum conservation and properties of 
C/ia/97- The tensor e^ Q/ g 7 has unit elements and is fully anti-symmetric: if 
£i234 = 1, other elements are obtained by changing the sign when adjacent 
elements are interchanged: e.g. €1243 = — 1. All elements with two identical 
indices are 0. 

To take account of the Breit-Wigner amplitude for the p and centrifugal 
barriers, should be multiplied by 

B 1 (k^)B 1 (k~) B 1 (k°-)B 1 (k+) B 1 (k-+)B 1 (k°) 

U (M 2 - s +0 - iMT) (M 2 - s - - iMT) (M 2 - s_ + - iMV) ' 1 ' 
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In the first term, B\(k ) is the L = 1 centrifugal barrier factor 



Bi{k~) = l/[k 2 _ + 0.03894/ J R 2 (/m)] 1 / 2 



(7) 



where k- is the momentum in GeV/c of the ir~ in the uj rest frame; R ~ 0.8 
fm. Also k +0 denotes the momentum of either ir in the p +0 rest system. The 
factor / w produces a mild modulation of the intensity of uj — > 7r + 7r~7r° decays 
over its Dalitz plot, but has no effect on angular correlations. 

In the uj rest frame, + + &o) 7 has only a time component, so W 7 ^ 
reduces to e^^fc+fc-M. The factor M for a; mass is not essential; remaining; 
factors reduce to the 3-dimensional vector product e^apk+k^L = k + Ak_; so W M 
is along the normal to the uj decay plane. It gives a weight factor \k + k-\f u . 
In the Y rest frame, from eqn. (5) does have a small time-component. 



form factors for the decay may have more effect than the variation of M over 
the resonance, so it can be dropped; BW(Y) is the Breit-Wigner amplitude 
for y(3941). Later expressions will omit the factor BW and centrifugal 
barrier factors, which are implied. In principle some centrifugal barrier factor 
is required for t = 1 in the production process. However, this form factor is 
likely to vary negligibly over the width of y(3941) and can be ignored. 

The full matrix element for B — > 1 ++ , 1 ++ — > [uj J/^]l=q is then the scalar 
product K.e A W = e.W A K. Here e has only x and y components in my 
axes, say (cos R, sin R). The intensity is given by the average over R, using 
< cos 2 R >=< sin 2 i? >= 1/2. In order to avoid this integration, a simple 
trick may be used in computer programmes. The e x and e y components may 
be replaced by 1 and i = Then intensities are obtained by taking the 

modulus squared of expressions like v a . 



For the production process, £ = 0, so the matrix element is simply the scalar 
product M = e a W a . One should form the 4-vector product, though both e± 



4.3 1++ [uJ/*] s=ljL=0 



The amplitude for this process in the Y rest frame is 



v, = e^ 4 e a W^M 4 BW(Y) = e^e a W^BW(Y); 



(8) 



4.4 ++ [uJ/y] s=0 ,L=o 
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and W4 are numerically small: 



P = ei W l + e 2 W 2 + e 3 W 3 - e 4 W A . (9) 

4.5 2++ -> [a;J/^] s=2 , L=0 

The fully relativistic expressions for tensors are given in Ref. [32]. Here, they 
are simplified to the rest frame of Y . For spin 2, one needs the tensor 

T aP = e a W p + e p W a - ^5 aP (e^), (10) 

where a, [3 and \i run from 1 to 3. This last term eliminates the scalar term 
which would otherwise be mixed into T a/ g; it makes the tensor traceless. 
The I = 2 operator for the production process is likewise constructed as 

r a p = K a Kp- l -5 a p{K 2 ). (11) 

As a reminder, K is the 3-momentum of the kaon in the rest frame of Y. 
The required matrix element for production and decay is then M = T a/3 T al3 . 
There is valuable information contained in the dependence on production 
angles through T a p. However, after integrating over all K directions, one 
finds the same decay angular distributions as given by Rosner's eqn. (17). 
This result arises from the fact that da/dVt = M*M contains a term r Q/3 r /3 ' Q! '. 
Integrating over angles, this term becomes 5 aa i5ppi for all amplitudes. 

4.6 Further amplitudes 

After these preliminaries, it is easy to write down remaining expressions. 
That for 1++ -> [u J/^] s =i, L =2 is 

M(l++ - [u;J/*] s=1 , i=2 ) = iTYV 3 , (12) 

where v 13 is given by eqn. (8) and r' is the L = 2 operator 

7 J afi = j a Jfi-\s afi (f l y t (is) 
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as a reminder, J is the 3-momentum of the J/^f in the rest frame of Y. Here 
the centrifugal barrier B 2 for L = 2 is needed: 

B 2 = 1/[(J 2 + Z)J 2 + Z 2 ] 1 / 2 , (14) 
Z = 0.11682/ R 2 (fm). 

Using both this and the L = amplitude eqn. (8) in the fit to F(3941) 
adds freedom which will need to be carefully controlled. The signature for 
the L = 2 amplitude is an interference with the S-wave proportional to J 2 . 
If the X(3872) and F(3941) have the same quantum numbers, they will have 
orthogonal wave functions, hence orthogonal amplitudes to a single channel, 
e.g. DD* or uJ/^. 

The amplitude for 1 ++ decaying to s = 2, L = 2 is less likely, since s = 2 
does not mix with s = 1, L = 0. It gives the same decay angular distribution 
after integrating over K and J, but may be ambiguous with s — 1, L — 1 for 
low statistics. It is the prototype for combining spins 2 and 2 to make J = 1: 

M(l++ - M/^] s=2iL=2 ) = (15) 

ee 2e a/37 T^V';. (16) 
For J PC< = 2 _+ , there are two amplitudes with s = 1 and 2: 

M(2- + - M/tt] a= i )L= i) = T a P[v a J p + v p J a -^5 a/3 (v.j)\, (17) 

M(2-+ - M/*] S=2 , L=1 ) = i^VZp (18) 

V^9 = e a ^rM 4 + e^ 4 T^rM 4 

= (e a ^ + 6^TZ)r. (19) 

This is the prototype for combining spins 2 and 1 to make 2. 
The amplitude for CT + is given by the simple contraction 

M(0- + -> [wJ/*] a= i )L= i) = J a t; a . (20) 

For J PC = 2++, there are three L = 2 amplitudes with s = 0, 1 and 2. 
The first with s = is given by the simple tensor contraction T a ^r' a ^. The 
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second for s = 1 is obtained from eqn. (19) by replacing J by v and T by r' . 
The third is 

M(2++ - [uJ/V] s=2 , L=2 ) = r^Y aP (21) 

Yap = TXd + T^-p^T^T'^). (22) 

There is little reason to expect these amplitudes to compete with L = 0. 
All four have the same dependence on the angle between lepton and W, and 
this will hopefully separate them from other spin-parities. The separation 
of possible s = 0, 1 and 2 amplitudes with L = 2 may be only of academic 
interest. 

The 3 — > [SJ/^]l=2 amplitude constructed from s — 1, L — 2 is 

2 

v itv\ = r V e A + T'n\e v + r'vxe^ - -{j 1 ^ a 5 vX + r' va 5 Xfl + r' Xa 5^)e a . (23) 

To form the full amplitude, this needs to be contracted with the operator for 
production with £ = 3: 

r, vX = K a K f3 K 1 - ^{K 2 ){K,5 vX + KJ^ + K x 5^). (24) 
4.7 Other channels 

Amplitudes for SJ/ty may be constructed along the same lines. For J PC = 
l + ~, one needs to combine e with J to make total spin 1, as in eqn. (8), 
replacing W by J and replacing p by the 5 amplitude. For J PC = 2 one 
combines s = 1 with L = 2, in analogy to eqn. (19), replacing J by e and T 
by r' . Amplitudes for 1 — > SJ/^ with L = and 2 are: 

M(l- -> [S-J/^^o) = K a e a (25) 
M(l~ - [SJ/*]l= 2 ) = K a r' a ^. (26) 

The decay B —> K + J/* is also given by eqn. (25). Production of ^(3770) 
with decay to [SJ/^/]l=2 is given by eqn. (26); for decays to the nn D-wave, 
M = K a T a pe^ . There is some evidence for ip{2S) decays to 07r° and ui] 
[33] . These may be fitted by replacing e by W for decays via the uj and by 
k(Ki) — k(K 2 ) for decays via the 0. 
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Amplitudes for DD* decays may be constructed in an analogous fashion. 
For C + 1, it is necessary to take the combination DD* + D*D; for C = — 1, 
the combination DD* — D*D in needed. The vector describing D* — > Dix is 
given by the 4- vector V 1 = (k D - K) - (M 2 D /M 2 D ^)(h D + k w ). For D* -> D7, 
the vector e of the J/\& is replaced by the corresponding vector for the photon. 

Amplitudes for decay to 7% are formed using the photon direction instead 
of the lepton from J/ty decays. The vector e so formed then needs to be 
combined with a vector constructed from xi decays or a tensor for \2 decays. 

A check on programming and amplitudes is that they should be orthogo- 
nal after integrating over all space. If formulae for higher spins are required, 
or for decays of B*, the methods described by Zou and Bugg [34] and Chung 
[35] are useful sources. 

5 Summary 

In Section 2, it has been shown that X(3872) decay via the tctt S-wave gives 
a poor fit to Belle data; J PC = 1 is not completely ruled out but would 
point to the vector glueball suggested by Seth [24]. Quantum numbers l + ~, 
2 and 3 are excluded by the poor fits with L — 1 and 2. 

Decays to pJ/^ with J PC = ++ , 1 ++ or 2 ++ are strong candidates; 
CT + and 2~ + are less likely alternatives. The latter could be eliminated if 
the 7T7T mass spectrum remains unchanged with a factor 4 higher statistics. 
They also have characteristic angular dependence discussed below. Higher 
statistics would also discriminate for or against 1 — > [SJ/^/]l =0 . 

Suggestions will be made here on how angular momentum analysis might 
proceed. Amplitudes with different J p are orthogonal, though experimental 
cuts might spoil this orthogonality to some extent. The X(3872) is narrow 
and presumably a single state. The dependence of amplitudes on five angles is 
distinctive and a separation between different J p with forseeable statistics of 
200 events looks realistic. For F(3941), more than one J p might contribute, 
making analysis more difficult. 

Amplitudes for F(3941) — > uJ/ty depend on the vector W a , eqn. (5), 
which lies normal to the decay plane of the 00. Amplitudes for AT(3872) — > 
pJ/^/ depend on the vector Q a = [k(ir + ) — k(n^)] a of eqn. (4). 

Firstly, the amplitude for J PC = 1 — > [SJ/^] L=0 is K.e. Since e lies 
normal to the direction of the lepton from J/ty decay, the amplitude is a 
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maximum with K normal to that direction. The intensity depends on sin 2 a, 
where a is the angle between kaon and lepton. Using the mass spectrum in 
addition, it should be possible to discriminate for or against this possibility. 

For pJ/ty and uJ/ty, decay angular distributions are distinctively differ- 
ent for J = ++ , f ++ and 2 ++ , as shown by Rosner [20]. Denoting by 
the angle between the lepton from J/ty decay and W (or Q), the decay an- 
gular distribution is proportional to sin 2 for ++ , (1 + cos 2 0) for 1 ++ and 
(7 — cos 2 0) for 2 ++ . The ++ amplitude is independent of the kaon direction, 
since there is no orbital angular momentum in the production process. 

For 1 ++ , the full matrix element is proportional to the scalar product 
e.W A K or e.Q A K for decays with L = 0. The amplitude is a maximum 
with the kaon direction perpendicular to W (or Q). There is then a sin 2 /? 
dependence on the angle (3 between the kaon direction and K or Q. Fur- 
thermore, W A K or Q A K will have a sine squared intensity variation with 
respect to the lepton direction. The overall angular dependence is much more 
distinctive than that on alone. If L — 2 decays of V(3941) also contribute, 
it is essential to constrain the analysis with the dependence of this amplitude 
on J 2 , eqn. (13); here J is the 3-momentum of the J/\I> in the rest frame of 
Y. The decay angular distribution is still proportional to (1 + cos 2 0). 

For J PC = 2 ++ , there are two distinctive features. The first is the decay 
angular distribution (7 — cos 2 0) , which is the same for all 2 + amplitudes. 
The second feature is the dependence on production through the tensor r a p of 
eqn. (11). For F(3941), the available momentum in the decay is 282 MeV/c 
at the peak. The suppression of the L = 2 amplitude by a conventional 
centrifugal barrier of radius 0.8 fm is then a factor 3 in amplitude, 9 in 
intensity. It is therefore likely that L = decays would dominate. The first 
step is to look for this amplitude r Q/3 T Q/3 where T is given by eqn. (10). The 
intensity has a fourth order dependence on the angle of the kaon with respect 
to the lepton from J/ty decay and a second order dependence on the angle 
between W and the lepton. This correlation is distinctive. If such angular 
dependence is observed, additional L = 2 amplitudes can be tried one by 
one. With modest statistics, it is not a good idea to fit using all four 2 + 
amplitudes, since experience elsewhere is that linear combinations of these 
four amplitudes can usually simulate lower spins 1 and 0. High statistics 
would be needed to separate the three L = 2 amplitudes with combined 
spins s = 0, 1 and 2 between u and J/ty and this may not be of great 
interest. 
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The J PC = h amplitude, eqn. (20), may be written as the scalar 
product e.J A W. The matrix element is a maximum when J and W are 
orthogonal and their vector product is normal to the lepton direction. The 
intensity has a sine squared dependence on the angle between J and W and 
on the angle between their vector product and the lepton. The amplitude 
is unique and angular correlations are very distinctive, so identifying this 
amplitude should be easy. 

For J PC = 2~ + there are two amplitudes, making analysis more difficult. 
As for 2 ++ , a characteristic signature is the dependence on T a/ 3, hence a fourth 
order dependence of intensity on the angles between the kaon and W, J or 
the lepton direction. The analysis is most easily done by fitting eqns. (17-19) 
to the data individually, then in combination. Again one should beware of 
two 2~ + amplitudes simulating lower spins. 

Now consider other decays. Information on the DD* channel is of primary 
importance. Formulae given here apply also to this channel and also decays 
to 7X and Sr) c , with simple substitutions of variables. If the X(3872) is a 
1 ++ bound state, sufficient statistics must reveal decays to DD* peaking in 
the first 10 MeV above threshold (in analogy with the NN 3 Si cross section, 
which is related to the existence of the deuteron). Evidence for 1 ++ D-wave 
decays would favour a bound-state; their absence would favour a cusp. The 
angular momentum analysis is most easily done by fitting S and D-wave 
1 ++ amplitudes directly to data. Decays through the D-wave would appear 
through interference with the S-wave and a distinctive dependence on J 2 , 
where J describes the direction of D or in the rest frame of DD*. This 
requires good mass resolution and quite high statistics. 

Charmonium radial excitations with J PC = ++ or 2 ++ are likely to have 
decays to DD and D S D S , so information on these channels is important. For 
spin 0, production and decay are isotropic. For spin 2, there is an intensity 
dependence (3 cos 2 7 — l) 2 , where 7 is the angle between D and the direction 
of the recoil kaon in the rest frame of Y . However, the weak process B — > 
K + 2 ++ is first forbidden, so the 2 ++ state may be suppressed there. 

Searches for open channels such as [St) c ]l=i are important. A knowledge 
of the rate of this process and that to DD* would allow a proper calculation 
of the magnitude to be expected from a cusp. 
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7 Appendix 



In the J/^f rest frame, the component of J/\& polarisation along the lep- 
ton axis is 0. Here, this polarisation is transformed first to Rosner's axes. 
Then the effect of the Lorentz transformation to the rest frame of F(3941) 
is evaluated. 

In my axes, e has components (cos/?, sin/?, 0), where R is unknown and 
must be averaged from to 2n. In Rosner's axes, e has components: 



cos 9 — sin 9 

1 
sin 9 cos 9 

I cos cos 9 cos R — sin sin R 
sin cos 9 cos R + cos sin R 
\ — sin 9 cos R 




(27) 
(28) 



Since e is a vector, the Lorentz transformation to the rest frame of y(3941) 
changes the Z-component to —7 sin 9 cos R and produces a time-component 
^7 sin 9 cos R. Rotating back through angles and 9 to axes in the Y rest 
frame parallel to the original z-axis, the result is a 4-vector 



( [1 + (7- l)sin 2 #]cosi? \ 
sini? 

sin 9 cos #(1 — 7) sini? 
/?7 sin 9 cos 9 cos R 



(29) 



For both F(3941) and X(3872), (3 is small, so the Lorentz transformation 
has only a small effect on e. 
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